De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

PeNAC67-PeKAN2-PeSCL23 and B-class MADS-box transcription factors synergistically regulate the specialization process from petal to lip in Phalaenopsis equestris.

The molecular basis of orchid flower development involves a specific regulatory program in which MADS-box transcription factors play a central role. The recent 'perianth code' model hypothesizes that two types of higher-order heterotetrameric complexes, namely SP complex and L complex, play pivotal roles in the orchid perianth organ formation. Therefore, we explored their roles and searched for other components of the regulatory network.Through the combined analysis for transposase-accessible chromatin with high-throughput sequencing and RNA sequencing of the lip-like petal and lip from Phalaenopsis equestris var.trilip, transcription factor-(TF) genes involved in lip development were revealed. PeNAC67 encoding a NAC-type TF and PeSCL23 encoding a GRAS-type TF were differentially expressed between the lip-like petal and the lip. PeNAC67 interacted with and stabilized PeMADS3, which positively regulated the development of lip-like petal to lip. PeSCL23 and PeNAC67 competitively bound with PeKAN2 and positively regulated the development of lip-like petal to petal by affecting the level of PeMADS3. PeKAN2 as an important TF that interacts with PeMADS3 and PeMADS9 can promote lip development. These results extend the 'perianth code' model and shed light on the complex regulation of orchid flower development.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives.

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Direct Bilirubin, but not Indirect Bilirubin, is Associated with Short-term Adverse Events in HFpEF.

OBJECTIVE: Abnormal live function tests have been identified as independent risk factors for ominous prognosis in patients with heart failure. However, most of the previous studies have failed to determine the contribution of direct bilirubin (DBIL) and indirect bilirubin (IBIL) separately. Hence, we aimed to explore whether DBIL or IBIL is correlated with the prognosis of heart failure with preserved ejection fraction (HFpEF).

Methods: A total of 19837 patients were hospitalized for HFpEF between January 2012 and January 2022 in Fuqing City Hospital affiliated with Fujian Medical University. The primary endpoint was in-hospital all-cause mortality. Secondary endpoints included in-hospital cardiovascular mortality and 30-day re-admission for heart failure.

Results: Univariable analysis indicated that patients with elevated DBIL or IBIL were exposed to a higher risk of mortality and re-admission. However, in multivariable models, both ln-transformed DBIL and TBIL, but not IBIL, were independent risk factors for in-hospital all-cause mortality [hazard ratio (HR)=1.796, 95% confidential interval (CI)=1.477-2.183, P＜0.001; HR=1.854, 95% CI=1.461-2.352, P＜0.001; HR=1.161, 95% CI=0.959-1.407, P=0.126] and in-hospital cardiovascular mortality (HR=1.831, 95% CI=1.345-2.492, P＜0.001; HR=1.899, 95% CI=1.300-2.773, P=0.001; HR=1.145, 95% CI=0.841-1.561, P=0.389). Only DBIL remained independently associated with 30-day readmission for heart failure (HR=1.361, 95% CI=1.036-1.787, P=0.027). Adding ln-transformed DBIL to model 1 increased its discriminatory capacity (C-statistic: 0.851 to 0.869, respectively), whereas adding ln-transformed IBIL yielded little increment (C-statistic: 0.851 to 0.852, respectively).

Conclusion: DBIL, but not IBIL, was associated with short-term ominous prognosis in patients with HFpEF. Hence, DBIL may be the superior predictor for prognosis in HFpEF.

AnimalMetaOmics: a multi-omics data resources for exploring animal microbial genomes and microbiomes.

The Animal Meta-omics landscape database (AnimalMetaOmics, https://yanglab.hzau.edu.cn/animalmetaomics#/) is a comprehensive and freely available resource that includes metagenomic, metatranscriptomic, and metaproteomic data from various non-human animal species and provides abundant information on animal microbiomes, including cluster analysis of microbial cognate genes, functional gene annotations, active microbiota composition, gene expression abundance, and microbial protein identification. In this work, 55 898 microbial genomes were annotated from 581 animal species, including 42 924 bacterial genomes, 12 336 virus genomes, 496 archaea genomes and 142 fungi genomes. Moreover, 321 metatranscriptomic datasets were analyzed from 31 animal species and 326 metaproteomic datasets from four animal species, as well as the pan-genomic dynamics and compositional characteristics of 679 bacterial species and 13 archaea species from animal hosts. Researchers can efficiently access and acquire the information of cross-host microbiota through a user-friendly interface, such as species, genomes, activity levels, expressed protein sequences and functions, and pan-genome composition. These valuable resources provide an important reference for better exploring the classification, functional diversity, biological process diversity and functional genes of animal microbiota.

SoyMD: a platform combining multi-omics data with various tools for soybean research and breeding.

Advanced multi-omics technologies offer much information that can uncover the regulatory mechanisms from genotype to phenotype. In soybean, numerous multi-omics databases have been published. Although they cover multiple omics, there are still limitations when it comes to the types and scales of omics datasets and analysis methods utilized. This study aims to address these limitations by collecting and integrating a comprehensive set of multi-omics datasets. This includes 38 genomes, transcriptomes from 435 tissue samples, 125 phenotypes from 6686 accessions, epigenome data involving histone modification, transcription factor binding, chromosomal accessibility and chromosomal interaction, as well as genetic variation data from 24 501 soybean accessions. Then, common analysis pipelines and statistical methods were applied to mine information from these multi-omics datasets, resulting in the successful establishment of a user-friendly multi-omics database called SoyMD (https://yanglab.hzau.edu.cn/SoyMD/#/). SoyMD provides researchers with efficient query options and analysis tools, allowing them to swiftly access relevant omics information and conduct comprehensive multi-omics data analyses. Another notable feature of SoyMD is its capability to facilitate the analysis of candidate genes, as demonstrated in the case study on seed oil content. This highlights the immense potential of SoyMD in soybean genetic breeding and functional genomics research.

Interaction between phenylpropane metabolism and oil accumulation in the developing seed of Brassica napus revealed by high temporal-resolution transcriptomes.

BACKGROUND: Brassica napus is an important oilseed crop providing high-quality vegetable oils for human consumption and non-food applications. However, the regulation between embryo and seed coat for the synthesis of oil and phenylpropanoid compounds remains largely unclear. RESULTS: Here, we analyzed the transcriptomes in developing seeds at 2-day intervals from 14 days after flowering (DAF) to 64 DAF. The 26 high-resolution time-course transcriptomes are clearly clustered into five distinct groups from stage I to stage V. A total of 2217 genes including 136 transcription factors, are specifically expressed in the seed and show high temporal specificity by being expressed only at certain stages of seed development. Furthermore, we analyzed the co-expression networks during seed development, which mainly included master regulatory transcription factors, lipid, and phenylpropane metabolism genes. The results show that the phenylpropane pathway is prominent during seed development, and the key enzymes in the phenylpropane metabolic pathway, including TT5, BAN, and the transporter TT19, were directly or indirectly related to many key enzymes and transcription factors involved in oil accumulation. We identified candidate genes that may regulate seed oil content based on the co-expression network analysis combined with correlation analysis of the gene expression with seed oil content and seed coat content. CONCLUSIONS: Overall, these results reveal the transcriptional regulation between lipid and phenylpropane accumulation during B. napus seed development. The established co-expression networks and predicted key factors provide important resources for future studies to reveal the genetic control of oil accumulation in B. napus seeds.

INSurVeyor: improving insertion calling from short read sequencing data.

Insertions are one of the major types of structural variations and are defined as the addition of 50 nucleotides or more into a DNA sequence. Several methods exist to detect insertions from next-generation sequencing short read data, but they generally have low sensitivity. Our contribution is two-fold. First, we introduce INSurVeyor, a fast, sensitive and precise method that detects insertions from next-generation sequencing paired-end data. Using publicly available benchmark datasets (both human and non-human), we show that INSurVeyor is not only more sensitive than any individual caller we tested, but also more sensitive than all of them combined. Furthermore, for most types of insertions, INSurVeyor is almost as sensitive as long reads callers. Second, we provide state-of-the-art catalogues of insertions for 1047 Arabidopsis Thaliana genomes from the 1001 Genomes Project and 3202 human genomes from the 1000 Genomes Project, both generated with INSurVeyor. We show that they are more complete and precise than existing resources, and important insertions are missed by existing methods.

Comparative transcriptome profiling reveals the multiple levels of crosstalk in phytohormone networks in Brassica napus.

Plant hormones are the intrinsic factors that control plant development. The integration of different phytohormone pathways in a complex network of synergistic, antagonistic and additive interactions has been elucidated in model plants. However, the systemic level of transcriptional responses to hormone crosstalk in Brassica napus is largely unknown. Here, we present an in-depth temporal-resolution study of the transcriptomes of the seven hormones in B. napus seedlings. Differentially expressed gene analysis revealed few common target genes that co-regulated (up- and down-regulated) by seven hormones; instead, different hormones appear to regulate distinct members of protein families. We then constructed the regulatory networks between the seven hormones side by side, which allowed us to identify key genes and transcription factors that regulate the hormone crosstalk in B. napus. Using this dataset, we uncovered a novel crosstalk between gibberellin and cytokinin in which cytokinin homeostasis was mediated by RGA-related CKXs expression. Moreover, the modulation of gibberellin metabolism by the identified key transcription factors was confirmed in B. napus. Furthermore, all data were available online from http://yanglab.hzau.edu.cn/BnTIR/hormone. Our study reveals an integrated hormone crosstalk network in Brassica napus, which also provides a versatile resource for future hormone studies in plant species.

BnIR: A multi-omics database with various tools for Brassica napus research and breeding.

In the post-genome-wide association study era, multi-omics techniques have shown great power and potential for candidate gene mining and functional genomics research. However, due to the lack of effective data integration and multi-omics analysis platforms, such techniques have not still been applied widely in rapeseed, an important oil crop worldwide. Here, we report a rapeseed multi-omics database (BnIR; http://yanglab.hzau.edu.cn/BnIR), which provides datasets of six omics including genomics, transcriptomics, variomics, epigenetics, phenomics, and metabolomics, as well as numerous "variation-gene expression-phenotype" associations by using multiple statistical methods. In addition, a series of multi-omics search and analysis tools are integrated to facilitate the browsing and application of these datasets. BnIR is the most comprehensive multi-omics database for rapeseed so far, and two case studies demonstrated its power to mine candidate genes associated with specific traits and analyze their potential regulatory mechanisms.

CottonMD: a multi-omics database for cotton biological study.

Cotton is an important economic crop, and many loci for important traits have been identified, but it remains challenging and time-consuming to identify candidate or causal genes/variants and clarify their roles in phenotype formation and regulation. Here, we first collected and integrated the multi-omics datasets including 25 genomes, transcriptomes in 76 tissue samples, epigenome data of five species and metabolome data of 768 metabolites from four tissues, and genetic variation, trait and transcriptome datasets from 4180 cotton accessions. Then, a cotton multi-omics database (CottonMD, http://yanglab.hzau.edu.cn/CottonMD/) was constructed. In CottonMD, multiple statistical methods were applied to identify the associations between variations and phenotypes, and many easy-to-use analysis tools were provided to help researchers quickly acquire the related omics information and perform multi-omics data analysis. Two case studies demonstrated the power of CottonMD for identifying and analyzing the candidate genes, as well as the great potential of integrating multi-omics data for cotton genetic breeding and functional genomics research.

High-throughput phenotyping-based quantitative trait loci mapping reveals the genetic architecture of the salt stress tolerance of Brassica napus.

Salt stress is a major limiting factor that severely affects the survival and growth of crops. It is important to understand the salt stress tolerance ability of Brassica napus and explore the underlying related genetic resources. We used a high-throughput phenotyping platform to quantify 2111 image-based traits (i-traits) of a natural population under three different salt stress conditions and an intervarietal substitution line (ISL) population under nine different stress conditions to monitor and evaluate the salt stress tolerance of B. napus over time. We finally identified 928 high-quality i-traits associated with the salt stress tolerance of B. napus. Moreover, we mapped the salt stress-related loci in the natural population via a genome-wide association study and performed a linkage analysis associated with the ISL population, respectively. These results revealed 234 candidate genes associated with salt stress response, and two novel candidate genes, BnCKX5 and BnERF3, were experimentally verified to regulate the salt stress tolerance of B. napus. This study demonstrates the feasibility of using high-throughput phenotyping-based quantitative trait loci mapping to accurately and comprehensively quantify i-traits associated with B. napus. The mapped loci could be used for genomics-assisted breeding to genetically improve the salt stress tolerance of B. napus.

Comparison of dynamic 3D chromatin architecture uncovers heterosis for leaf size in Brassica napus.

INTRODUCTION: Heterosis is the major event driving plant development and promoting crop breeding, but the molecular bases for this phenomenon remain elusive. OBJECTIVES: We aim to explore the effect of three-dimensional (3D) chromatin architecture on the underlying mechanism of heterosis. METHODS: Here, we constructed the North Carolina II (NC-II) population to select superior and inferior heterosis sets by comparing mid-parent heterosis (MPH) in Brassica napus. To decipher the impact of 3D chromatin architecture on the underlying mechanism of heterosis, we combined genetics, transcriptomics and 3D genomics approaches. RESULTS: We suggest that F1 hybrids with superior heterosis tend to contain more transcriptionally active A compartments compared with F1 hybrids with inferior heterosis, and approximately 19-21% compartment significantly altered in the F1 hybrids relative to the parental lines. Further analyses show that chromatin compartments correlate with genetic variance among parents, which may form the basis for differentially active chromatin compartments. Having more A compartments in F1 hybrids confers a more accessible chromatin circumstance, which promotes a higher proportion of highly expressed ELD (expression level dominance) genes in superior heterosis F1 hybrids (46-64%) compared with inferior heterosis F1 hybrids (22-31%). Moreover, genes related to hormones which affect plant growth, are more up-regulated with changes of 3D genome architecture, and we validate that increased hormone content contributes to cell proliferation and expansion by influencing the key genes of cell cycle thereby promoting leaf size. CONCLUSION: Dynamic 3D chromatin architecture correlates with genetic variance among parents and contributes to heterosis in Brassica napus.

Asymmetric gene expression in grain development of reciprocal crosses between tetraploid and hexaploid wheats.

Production of viable progeny from interploid crosses requires precise regulation of gene expression from maternal and paternal chromosomes, yet the transcripts contributed to hybrid seeds from polyploid parent species have rarely been explored. To investigate the genome-wide maternal and paternal contributions to polyploid grain development, we analyzed the transcriptomes of developing embryos, from zygote to maturity, alongside endosperm in two stages of development, using reciprocal crosses between tetraploid and hexaploid wheats. Reciprocal crosses between species with varied levels of ploidy displayed broad impacts on gene expression, including shifts in alternative splicing events in select crosses, as illustrated by active splicing events, enhanced protein synthesis and chromatin remodeling. Homoeologous gene expression was repressed on the univalent D genome in pentaploids, but this suppression was attenuated in crosses with a higher ploidy maternal parent. Imprinted genes were identified in endosperm and early embryo tissues, supporting predominant maternal effects on early embryogenesis. By systematically investigating the complex transcriptional networks in reciprocal-cross hybrids, this study presents a framework for understanding the genomic incompatibility and transcriptome shock that results from interspecific hybridization and uncovers the transcriptional impacts on hybrid seeds created from agriculturally-relevant polyploid species.

The genome of Dioscorea zingiberensis sheds light on the biosynthesis, origin and evolution of the medicinally important diosgenin saponins.

Diosgenin saponins isolated from Dioscorea species such as D. zingiberensis exhibit a broad spectrum of pharmacological activities. Diosgenin, the aglycone of diosgenin saponins, is an important starting material for the production of steroidal drugs. However, how plants produce diosgenin saponins and the origin and evolution of the diosgenin saponin biosynthetic pathway remain a mystery. Here we report a high-quality, 629-Mb genome of D. zingiberensis anchored on 10 chromosomes with 30 322 protein-coding genes. We reveal that diosgenin is synthesized in leaves ('source'), then converted into diosgenin saponins, and finally transported to rhizomes ('sink') for storage in plants. By evaluating the distribution and evolutionary patterns of diosgenin saponins in Dioscorea species, we find that diosgenin saponin-containing may be an ancestral trait in Dioscorea and is selectively retained. The results of comparative genomic analysis indicate that tandem duplication coupled with a whole-genome duplication event provided key evolutionary resources for the diosgenin saponin biosynthetic pathway in the D. zingiberensis genome. Furthermore, comparative transcriptome and metabolite analysis among 13 Dioscorea species suggests that specific gene expression patterns of pathway genes promote the differential evolution of the diosgenin saponin biosynthetic pathway in Dioscorea species. Our study provides important insights and valuable resources for further understanding the biosynthesis, evolution, and utilization of plant specialized metabolites such as diosgenin saponins.

Genome-wide association study and transcriptome analysis reveal key genes controlling fruit branch angle in cotton.

Fruit branch angle (FBA), a pivotal component of cotton plant architecture, is vital for field and mechanical harvesting. However, the molecular mechanism of FBA formation is poorly understood in cotton. To uncover the genetic basis for FBA formation in cotton, we performed a genome-wide association study (GWAS) of 163 cotton accessions with re-sequencing data. A total of 55 SNPs and 18 candidate genes were significantly associated with FBA trait. By combining GWAS and transcriptome analysis, four genes underlying FBA were identified. An FBA-associated candidate gene Ghi_A09G08736, which is homologous to SAUR46 in Arabidopsis thaliana, was detected in our study. In addition, transcriptomic evidence was provided to show that gravity and light were implicated in the FBA formation. This study provides new insights into the genetic architecture of FBA that informs architecture breeding in cotton.